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Abstract 



CZ5 

k^ V The present work investigates the direct interaction of sheared mean 

, &-f - flow with zonal flows (ZF) and the effect of parallel ion motion on ZF 

generation in ion-temperature-gradient (ITG) background turbulence. 

> , 

I/"") \ An analytical model for the direct interaction of sheared mean flows 

t-H ! 

■ with zonal flows is constructed. The model used for the toroidal ITG 

(N . 

driven mode is based on the equations for ion continuity, ion tempera- 

o . 

ON , ture and parallel ion motion whereas the ZF evolution is described by 



the vorticity equation. The behavior of the ZF growth rate and real 



^ ■ frequency is examined for typical tokamak parameters. It is shown 

that in general the zonal flow growth rate is suppressed by the pres- 
ence of a sheared mean flow. In addition, with parallel ion motion 
effects the ZFs become more oscillatory for increasing rji(= L n /LTi) 
value. 
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I Introduction 



The study of plasma flows for regulating the turbulence and anomalous trans- 
port has in recent years attracted strong interest [TJ. The regulation of drift 
wave turbulence by sheared flows is due to shearing of the turbulent eddies 
and thereby reducing the spatial scales of the eddies. Of particular inter- 
est among the plasma flows are the so called zonal flows [2] which are self 
generated from the background turbulence via Reynolds stress. 

Unlike the sheared mean flows which have equilibrium scale size the zonal 
flows are random E x B flows mainly in the poloidal direction with low 
frequency and thus are almost stationary compared to the time scale of the 
background turbulence. Since the zonal flows are quasistationary compared 
to the background turbulence they may keep decorrelating turbulent eddies 
for a relatively long time and thereby effectively suppressing the turbulent 
transport. 

The drift waves that are driven by gradients in the plasma density, tem- 
perature, magnetic field etc are responsible for causing turbulence and anoma- 
lous transport. Candidates of such drift waves are the Ion- Temperature- 
Gradient Mode (ITG) and the collisionless Trapped Electron Mode (TEM). 

In the present paper the direct effects of a sheared mean flow on the 
zonal flow instability are investigated analytically, however, the effects on the 
background turbulence itself is not considered in this work. Moreover, the 
effect of the parallel ion motion on the zonal flow growth rate and frequency 
is studied. 

The zonal flow generation from nonlinear interactions among drift waves 
has been extensively investigated both analytically [3]- [15] and in computer 
simulations using gyrokinetic [TTJ- [TH] and advanced fluid models [20]- [23J . 



2 



Furthermore, the effect of micro-scale electron temperature gradient (ETG) 
driven turbulence driven zonal flow on semi-macro scale ITG turbulence was 
studied by a fluid simulation [24J. However, the direct interaction of zonal 
flow with a sheared mean flow and the interaction with the GAM have mostly 
been overlooked. There have been some previous studies on the interaction 
of zonal flows and mean flows using simple drift wave models [T3j and using 
the coherent mode coupling method [16] . The present work extends the drift 
wave model to an advanced fluid model for the toroidal ITG mode turbulence 
while using the wave kinetic approach. The purpose of this study is to obtain 
a qualitative estimate of the zonal flow growth rate and real frequency, and 
their parametric dependence of the plasma under the influence of a sheared 
mean flow and parallel ion motion. The effects of sheared mean flows and 
parallel ion motion on zonal flows is not well known as well as how it is 
influencing the turbulence and in the end how this effect may change the 
overall transport. 

However the full non-linear effects of the Geodesic Acoustic Mode (GAM) [25] 
is out of the scope of the present study. The GAM is a perturbation where 
the electrostatic potential (m = n = mode) is coupled to the sideband 
density perturbations (m = l,n = mode) by toroidal effects. The GAM 
interacts with the zonal flow and act as a source or sink for the poloidal 
flow leading to an oscillatory nature of the zonal flow. The complete role of 
GAMs are still not well understood [26]- [29]. For these reasons it is of great 
importance to study the four component system of drift waves, sheared mean 
flow, GAM and zonal flow. 

The methodology of the analytical model is based on the wave kinetic 
approach [5], [13]- [15]. An algebraic equation which describes the zonal 
flow growth rate and real frequency including the effects of a sheared mean 
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flow in the presence of maternal ITG turbulence including the effects of 
parallel ion motion is derived and solved numerically. An advanced fluid 
model including the ion continuity, ion temperature equations and the ion 
momentum equation is used for the background ITG turbulence [5U]. The 
generation of zonal flows is described by the vorticity equation and the time 
evolution of the ITG turbulence in the presence of the slowly growing zonal 
flow is described by a wave kinetic equation. 

It was found the generation of zonal flows was in general suppressed by 
a sheared mean flow. By introducing collisional damping a modest suppres- 
sion of the zonal flow growth rate is found. In addition, it is found that 
the interaction with parallel iom momentum may reduce the ZF generation 
significantly. 

The paper is organized as follows. In Section II the analytical model for 
the zonal flows generated from toroidal ITG modes including the effects of a 
sheared mean flow is reviewed. The analytical model is extended to include 
parallel ion motion in Section III. Section IV is dedicated to the results and 
a discussion thereof. Finally there is a summary in section V. 

II Analytical model for interaction of mean 
flows and zonal flows 

In this section the model excluding parallel ion motion is introduced. The 
description used for toroidal ITG driven modes consists of the ion continuity 
and ion temperature equations. For simplicity, effects of electron trapping 
and finite beta are neglected in this work. Magnetic shear can, however, 
modify the non-linear upshift as found in Ref. [TU]- [H] and is accordingly 
incorporated in the model including parallel ion motion. In this section 
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the model of how to construct the interaction of a sheared flow and the 
zonal flows generated from ITG driven turbulence and the derivation of the 
dispersion relation for zonal flows are summarized. The method has been 
described in detail in Refs. [5], [T3]- [T5] (and References therein) and only a 
brief summary is given here. The analysis of Ref. [S] is closely followed and 
in the present work extended to be valid for an advanced fluid model. An 
alternative statistical approach, resulting in a modified wave kinetic equation, 
is presented in Ref. [12] which also contains an extensive discussion of and 
comparison with the approach used here. In describing the large scale plasma 
flow dynamics it is assumed that there is a sufficient spectral gap between the 
small scale fluctuations and the large scale flow. The electrostatic potential 
(0 = eip/T e ) is represented as a sum of fluctuating and mean quantities 

0(X,x,T,t) = $(X,T) + 0(f,t) (1) 

where (p(x, t) is the fluctuating potential varying on the turbulent scales x, y, t 
and $(X,T) is the zonal flow potential varying on the slow scale X,T (the 
zonal flow potential is independent on Y). The coordinates (XjTjj (x, t) 
are the spatial and time coordinates for the mean flows and small scale fluc- 
tuations, respectively. 

The wave kinetic equation see Refs. [5], [31] - [32] for the generalized 
wave action N k = A ^ fc 7 ^ 1 4>k \ 2 in the presence of a sheared mean plasma 
flow perturbing the other mean flow (the zonal flow in this case) due to the 
interaction between mean flow and small scale fluctuations is 

d M( 8 , t ^ dN k (x,t) d (t ^dN k (x,t) 



X 



= j k N k (x,t) - AujN k (x,t) 2 (2) 

Here the spectral difference (q y << k y , q y ,k y is zonal flow and drift wave 
wave numbers respectively) is used in solving the wave kinetic equation and 



only the x direction is considered where similar spectrum of the background 

turbulence and zonal flow is expected. Where vq is the zonal flow part of the 

E x B velocity, and in the relation between small scale turbulence and the 

generalized wave action density we have 
k ( At \ 

A fc = "7T ( 1 — e n9 + 1T^n9 (3) 



Ik = k y \Je n g (Vi-Vith)- (4) 

In the expression for the rjuh the FLR effects are neglected, 

21 1 ( 1 10 \ /.X 

^"3-27 + 4^ + e ^U + 97j- (5) 

Here, and in the forthcoming equations r = Tj/T e , v+ = p s c s y/L n , p s = c s /fl ci 



where c s = yT e /rrii, Vt ci = eB/rriiC. We also define Lf = — (dlnf /dr)^ 1 , 
r/i = Ln/L/Ti, e n = 2L n /R where R is the major radius and oii = r (1 + rji). 
The perturbed variables are normalized with the additional definitions n = 
(L n /p s )5n/n , (f> = (L n / p s e)5(f>/T e , % = (L n / p s )5Ti/T i0 as the normalized 
ion particle density, the electrostatic potential and the ion temperature, re- 
spectively. The perpendicular length scale and time are normalized to p s 
and L n /c s , respectively. The geometrical quantities are calculated in the 
strong ballooning limit (8 = 0, g (9 — 0, k) = 1/k (Ref. [33]) where g (6) 
is defined by uJn{0) = u^e n g(9)). In this analysis it is assumed that the 
RHS is approximately zero (stationary turbulence). The role of non-linear 
interactions among the ITG fluctuations (here represented by a non-linear 
frequency shift Au) is to balance linear growth rate. In the case when 
7fcA^fc(x,t) — AuNklxjt) 2 = 0, the expansion of the wave kinetic equation 
is made under the assumption of small deviations from the equilibrium spec- 
trum function; Nj. = N® + Nk where Nk evolves at the zonal flow time and 
space scale (Q,q x ,q y = 0), as 

-qT- + iqxV gx N k - k y (Ve) -qJ-~ + l k N k = iqx.k y V E -^- (6) 



In this last expression the third term on the left hand side shows the explicit 
interaction of a mean shear flow (V E ) and Nk- Here Ve = iq x ^- This 
equation is solved for Nk assuming that (V E ) 2 is small, by integrating by parts 
and obtaining an expansion in (V E ) 2 and introducing a total time derivative 
D t = j^ — ky (Vg) It is interesting to note that on the total time scale the 
shearing effect is explicit as D t k x = —k y (V E ). The solution can be written 
as 

N k = / dt'e-^-^-^Jf dtv "Hq x k y V E ^- 

J to OK x 

f)N° 

= - iq lk y ^R (k y ,q x ,n,(V E ))^. (7) 

The evolution equations for the zonal flows is obtained after averaging the 
ion-continuity equation over the magnetic flux surface and over fast scales 
employing the quasineutrality and including a damping term [31]. The 
evolution equation is obtained as 

|V.. - = (1 + r) VI + rVl » f tt ) (8) 

where it is assumed that only the small scale self interactions are the im- 
portant interactions in the RHS [35J. Using typical tokamak parameters 
{Ti = T e = lOkev, m = n e = KPm" 3 , r = lm, R = 3m) \i = Q.7%v iiy f(r/R) 
and Vu = 10~ 12 rii/T^ 2 and fa is the ion-ion collision frequency, Tj is the ion 
temperature in electron volts. Using typical tokamak parameters it is found 
that /i ~ 50. Expressing the Reynolds stress terms in Eq. 8 in Nk we obtain 

{-in - wl) $ = (1 + r + tS) J d 2 kk x k y \4> k \ 2 (9) 

where 5 is a k independent factor 

A k k v ( 2 , . \ , . 

J = Alt%("'-3 (1 + T) ^)- (10) 
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Utilize equations 7, 9 and 10 gives, 

A 2 -I- 'v 2 r f)M° 
(-<n - *£) = -<£ (1 + r + r*) / d 2 kk 2 k x ^-R Q . (11) 

Here the response function (i? ) is, considering only the first two even terms 
in the expansion 

Ro = c (12) 

7 fc - - v gx q x ) 

R = R + R 2 D t R D t . (13) 



In contrast with Ref. [14] it is not assumed that the short scale turbulence 
is close to marginal state (or stationary state, 7& is small). Integrating by 
parts in k x and assuming a monochromatic wave packet iV° = Nq5 (k — ko) 
gives 

(ft + iwl) = -iq 2 x (1 + r + t5) ^-^ik 2 y Ro(k y , q X) ft, (V E ))N (14) 

Ik 

Here, ft is the zonal flow growth rate and real frequency, q x is the zonal flow 
wave number and k y is the wave number for the ITG mode. The real part of 
the sheared mean flow dependent term can now be written in the k± << 1 
limit as found in Ref. [T3j . 

7fe 



Re(Rc 



7i + (ft - qxv gx ) 2 

iw ^V-^ ) 5 - (15) 



This result is valid in the weak shear limit 'Jzf > (Ye)- m later numerical cal- 
culations the full expression is retained. In expressing the zonal flow growth 
in dimensional form making use of the relation (A| + 7^) / (47^) N = |0| 2 it 
is assumed that the mode coupling saturation level is reached [36] 

7 1 
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When calculating the group velocities the FLR effects in linear ITG mode 
physics is important and given by the real frequency and growth rate as 
follows from Ref Hi 



u r = , 1 fcy , 2 , fl - ( 1 + tn9 - k\ (a, + %re n g ) ) (17) 



2 (1 + k\) \ V 3 J M x V 1 3 

k I 

7 = l + k\ v T€n9 ( 18 ) 

where u = oj r + i-y. The group velocities (v g j = du r /dkj) are in the long 

wavelength limit (k^_ << 1) given by, 



v 



(jX 



5r 



k x k y [1 + (1 +T]i)r — [1 + — ) e n g) (19) 



«» = + (20) 



III Analytical model including parallel ion mo- 
tion 

It is known that parallel ion motion effects on the background turbulence 
growth rate is only slightly modified, whereas, there is often a significant 
effect on the real frequency (significant increase in \u r \). The zonal flow 
dispersion equation is explicitly dependent of the background real frequency, 
group velocity and growth rate and it is now expected that there is significant 
effect on the zonal flow generation dependent on the parallel ion motion. 
Accordingly, the previous model for the generation of zonal flows from ITG 
background turbulence is extended to include the equation of motion for the 
ions. The model for the drift waves consists of the following equations: 
ion continuity equation 

dn ( d d \ o 7 36 d /- / ~\\ dv,\\ 



dt \dt dy J dy dy ^ ^ ' ' dz 

, n) + U, Vi0l + r U, Vi (n + Ti)] (21) 



ion energy equation 

dfi 5 dfi ( 2\d<p 2dn r(ml 2 r , . 



parallel ion momentum equation 



du 



dt \I + T ^(" + ^) 



d 



(23) 



Here [A, B] = dA/dxdB/dy — dA/dydB jdx is the Poisson bracket. The 
quantities are normalized in the same fashion as above with = [L n j p s )vi\\/c s 
The electrons are assumed to be Boltzmann distributed. Note that, for the 
zonal flows k\\ = and Eq. 23 is identically zero and Boltzmann distributed 
electrons cannot be used, instead the same model as earlier is employed c.f 
Eq. 8. The dispersion relation for the ITG mode resulting from Eqs 21-23 
is then 

5r 



1 + kl ( 14 



5r 

y 



1 - e n 1 + 



+ 



u 2 - 



-<w( i+ t); 



a r ) — k ,tT 



00 k,, 



e n I r - a r + — k ± (1 + 77.) J + % , , 



fej = 



(24) 



Here a r = ^, r = T e /Tj and 



r = r (rii - - J + —e n (1 + r) 



(25) 



Here, the solution have been found using an approximate eigen-mode function 
in the form of the lowest order Hermite polynomial (n — 0). 



5(f) oc e 



-z 2 /2a 2 



(26) 



where 



k±\s\qu 



(27) 
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Here, Re(cr~ 2 ) > has been imposed as a causality requirement. Now, we 
proceed as in Ref. [H] and obtains a modified 3rd order dispersion relation 
for the zonal flow growth rate and real frequency. 



where u r , 7^ and the numerical derivative |^ are found from the dispersion 
relation (Eq. 24), Rq is the same factor as before and where the factor S now 
is modified to, 



In present model extended for including parallel momentum and mode struc- 
ture there are some limitations, it is only a change in the linear drift wave 
physics and the possible non-linear coupling term in Eq. 28 is not included 
and analogously the generalized action density derived in [T3] for the back- 
ground turbulence is only approximately correct in the long wavelength limit. 

IV Results and discussion 

The model in Ref. [13] is expanded by including the effects of interaction 
with mean flows and the effects of parallel ion motion for the zonal flow 
generation. An analytical dispersion relation is derived according to the 
wave kinetic equation approach that takes into account the effect of a shear 
flow interacting with the generation of zonal flows. The plasma parameter 
dependence on the zonal flow growth rate and real frequency is explored 
for comparison purposes with the explicit dependence of mean flow shear 
{(Ye) 0(1 Q}^f) gi ven special attention. Here q$ is the wavenumber and $/ is 
the electrostatic potential of the sheared mean flow. First the results from 
the analytical model derived from the wave kinetic approach is presented 




(28) 




(29) 
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Figure 1: the zonal flow growth rate and frequency as a function of rji with 
collisional damping (/i) as a parameter. The other parameters are as q x = 
0.3 — k x — k y , e n = 1.0 and r = 1, qf — and |$/| = 0. The results are 
shown for jj, = (plus), ;U = 10 (boxes) and /i = 50 (rings). 

including the effect of a sheared mean flow (solution to the dispersion relation 
Eq. 12). 

In Figure 1 the zonal flow growth rate and frequency as a function of rji 
with the collisional damping as a parameter is displayed. The other parame- 
ters are q x = 0.3 = k x = k y , r = 1 and e n = 1 with qf = and \&f\ = 0. The 
results are shown for ji = (plus), ji = 10 (boxes) and /i = 50 (rings). In the 
figure the growth rates are positive while the real frequencies are negative for 
the zonal flow. In the case of no damping fi = 0, almost no explicit r\i depen- 
dency on the zonal flow growth rate and real frequency appears, whereas, for 
increasing damping (/i) modest effects on the growth rate and real frequency 
of the zonal flow f2 are exhibited. 

Next in Figure 2 and 3, the zonal flow growth rate (Figure 2) and real 
frequency (Figure 3) as a function of rji with |$/| as a parameter is shown. 
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The other parameters are as in Figure 1 with q$ = 0.1. The results are 
shown for |$/| = (plus), |$/| = 100 (boxes) and |$/| = 500 (asterisk). 
The effect of a mean flow on zonal flow generation is qualitatively similar to 
that of damping, however, in the mean flow case the growth rates and real 
frequencies approach the case with no flow |$j| = for large values of r\i. In 
addition, unlike the damping, the sheared mean flow is also responsible for 
a significant suppression of the zonal flow growth rate for small ^-values to 
a certain minimum level determined by the background turbulence. If the 
background is considered to be close to the marginal state a total suppression 
of the zonal flow growth rate is found. Note that the effect of mean flow 
is included in the product of qj and |$/|, namely the shearing rate. The 
resulting zonal flow growth rate reflects that the solution to the dispersion 
relation consists of two branches where one branch is not depending on the 
mean flow damping and another branch that is significantly damped close to 
marginal stability. The behavior of the zonal flow growth rate is indicated in 
Eq. 15 where the mean flow damping is explicitly found and that it is reduced 
with increasing r^- values (i.e. increasing ITG growth rate). If the sheared 
mean flow is allowed to interact with the background and act as E x B shear 
the linear ITG growth rate will be reduced resulting in a reduction of the 
zonal flow growth rate. Moreover, in Eq. 15 it is indicated that a reduction 
in the linear growth rate will also result in a stronger direct effect of a sheared 
mean flow damping on the zonal flows. 

In Figure 4 and 5 the zonal flow growth rate (Figure 4) and real frequency 
(Figure 5) (normalized to the linear ITG growth rate) as a function of e n is 
shown with mean flow shear as a parameter is displayed. The other param- 
eters are as in Figure 1 with r/i = 4 and qf — 0.1. The results are shown for 
|$/| = (asterisk), |$/| = 100 (diamonds) and |$/| = 500 (stars). Note, 
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Figure 2: (2 and 3) The zonal flow growth rate (Figure 2) and frequency 
(Figure 3) as a function of rji with (?/|$/| as a parameter. The other param- 
eters are as q x = 0.3 = k x = k y , e n = 1.0 and r = 1. The results are shown 
for <?/|<I>/ 1 = (plus), qf\®f\ = 100 (boxes) and = 500 (asterisk). 
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Figure 3: (2 and 3) The zonal flow growth rate (Figure 2) and frequency 
(Figure 3) as a function of rji with g/|3>/| as a parameter. The other param- 
eters are as q x = 0.3 = k x = k y , e n = 1.0 and r = 1. The results are shown 
for = (plus), 1 = 100 (boxes) and ?/|3>/| = 500 (asterisk). 



14 



0.5 



2 



Figure 4: (4 and 5) The zonal flow growth rate (Figure 4) and frequency 
(Figure 5) (normalized to the linear ITG growth rate) as a function of e n 
is shown with mean flow shear as a parameter. The other parameters are 
q x = 0.3 = k x = k y , r = 1 and r)i = 3. The results are displayed for 
<Z/|$/| = (asterisk), = 100 (diamonds) and g/l^/l = 500 (stars). 

that in the present normalization the expansion term (k y qj\§f\ |) is still small 
compared to the growth rate in the weak flow case while it is comparable to 
the growth rate in the strong flow case. The results show a suppression of 
zonal flow growth with |$/|. The reason for this is in the variation of the 
radial group velocity (v gx ) in the total time frame, as in the decorrelation of 
drift wave propagation by a sheared flow. This will weaken the modulation 
response of the drift wave response. The stabilization of zonal flow growth by 
the sheared flow results in a significant increase of the real frequency (|f2 r |). 
The effects resulting from from this approach is in qualitative agreement 
with with results reported earlier using other drift wave models [13J and us- 
ing the coherent mode coupling method in ETG background turbulence for 
generating the zonal flow [16J. 
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Figure 5: (4 and 5) The zonal flow growth rate (Figure 4) and frequency 
(Figure 5) (normalized to the linear ITG growth rate) as a function of e n 
is shown with mean flow shear as a parameter. The other parameters are 
q x = 0.3 = k x = k y , r = 1 and rji = 3. The results are displayed for 
qf\$f \ = (asterisk), = 100 (diamonds) and = 500 (stars). 
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Next, the results using the extended model including parallel ion momen- 
tum is treated (solutions to the dispersion relation Eq. 28). 

In Figure 6, the ratio of the growth rate and real frequency is shown as a 
function of rji with safety factor (q) as a parameter. In the present case the 
safety factor is varied from q = 2 (asterisk curve) to q = 8 (plus curve). The 
other parameters are as in Figure 1. It is found that the ratio of the growth 
rate and the frequency is decreasing with increasing rji (remember that L n 
is considered to be fixed). This is suggestive of a transition from a state of 
stable zonal flows to a state with oscillating zonal flows. It is indicated that 
in the region close to the linear ITG threshold the zonal flows are stationary 
and may have a significant stabilizing effect on the background turbulence 
whereas at higher rji the zonal flows becomes oscillatory. In general, in a 
comparison of the previous model and the present model the change in zonal 
flow generation inherently comes from the fact that \u r \ is increased in the 
system including parallel ion momentum, whereas the ITG growth rates are 
rather similar in both systems. 



V Summary 

This work focus mainly on analytical estimations of zonal flow growth rate 
and real frequency under the influence of sheared mean flows and parallel 
ion motion. An analytical model for the generation of zonal flows by ion- 
temperature-gradient background turbulence in the presence of a sheared 
mean flow and parallel ion motion are derived using the wave kinetic ap- 
proach. The model consists of the ion continuity and ion temperature equa- 
tions, in addition, the parallel momentum equation for the ions are included. 
The zonal flow evolution is described by the vorticity equation including a 



17 



Figure 6: The ratio of the growth rate and frequency 'Jzf/^zf as a function 
of r)i with safety factor (q) as a parameter. The other parameters are as 
q x = 0.3 = k x = k y , e n = 1.0, r = 1 and = 0. The results are 

displayed for q = 2 (asterisk curve) to q = 8 (plus curve). 

collisional damping on the zonal flow generation. The zonal flow growth rates 
and real frequencies are scanned for a wide range of plasma parameters. 

It was found the general level of zonal flow was suppressed by a sheared 
background flow and also that for obtaining realistic results retaining colli- 
sions into the analytical model seems to be very important. 

The results show a suppression of zonal flow growth with | $/ 1 . This is in 
agreement with previous research of zonal flow generation using another drift 
wave model [13] and the coherent mode coupling method in ETG drift wave 
turbulence [16J. The reason for this is in the variation of the radial group 
velocity (v gx ) in the total time frame, as in the decorrelation of drift wave 
propagation by a sheared flow. This will weaken the modulation response 
of the drift wave response. The stabilization of zonal flow growth by the 
sheared flow results in a significant increase of the real frequency (|n r |). 
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By introducing collisional damping a suppression of the zonal flow growth 
rate is found. In the case of no damping fi = 0, no explicit rji dependency 
on the zonal flow growth rate and real frequency appears, whereas, a modest 
effect of damping is exhibited. 

In addition, it is found that the parallel ion motion may reduce the ZF 
generation significantly. This is suggestive of a transition from a state of 
stable zonal flows to a state with oscillating zonal flows. It is indicated that 
in the region close to the linear ITG threshold the zonal flows are stable 
and may have a significant stabilizing effect on the background turbulence 
whereas at higher rji the zonal flows becomes oscillatory. In general, in a 
comparison of the previous model and the present model the change in zonal 
flow generation inherently comes from the fact that \uj r \ is increased in the 
system including parallel ion momentum, whereas the ITG growth rates are 
rather similar in both systems. 
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